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Abstract 

The Schrodinger equation in D-dimensions for the Manning-Rosen poten- 
tial with the centrifugal term is solved approximately to obtain bound states 
eigensolutions (eigenvalues and eigenfunctions) . The Nikiforov-Uvarov (NU) 
method is used in the calculations. We present numerical calculations of en- 
ergy eigenvalues to two- and four-dimensional systems for arbitrary quantum 
numbers n and / with three different values of the potential parameter a. It 
is shown that because of the inter dimensional degeneracy of eigenvalues, we 
can also reproduce eigenvalues of a upper/lower dimensional sytem from the 
well-known eigenvalues of a lower / upper dimensional system by means of the 
transformation (n,l,D) (n, Z it 1,Z) =p 2).. This solution reduces to the 
Hulthen potential case. 
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I. INTRODUCTION 

One of the important tasks of quantum mechanics is to find exact solutions of the wave 
equations (nonrelativistic and relativistic) for certain potentials of physical interest since 
they contain all the necessary information regarding the quantum system under consider- 
ation. It is well known that the exact solutions of these wave equations are only possible 
in a few simple cases such as the Coulomb, the harmonic oscillator, pseudoharmonic and 
Mie-type potentials [1-8]. For an arbitrary Z-state, most quantum systems could be only 
treated by approximation methods. For the rotating Morse potential some semiclassical 
and/or numerical solutions have been obtained by using Pekeris approximation [9-13]. In 
recent years, many authors have studied the nonrelativistic and relativistic wave equations 
with certain potentials for the s- and /-cases. The exact and approximate solutions of these 
models have been obtained analytically [10-14]. 

Many exponential-type potentials have been solved like the Morse potential [12,13,15], 
the Hulthen potential [16-19], the Poschl-Teller [20], the Woods-Saxon potential [21-23], the 
Kratzer-type potentials [12,14,24-27], the Rosen- Morse-type potentials [28,29], the Manning- 
Rosen potential [29-33] and other multiparameter exponential-type potentials [34,35]. Var- 
ious methods are used to obtain the exact solutions of the wave equations for this type of 
exponential potentials (cf. [36] and the references therein). 

Recently, the NU method [37] has shown its power in calculating the exact energy levels 
of all bound states for some solvable quantum systems. In this work, we attempt to apply 
this method to study another exponential-type potential proposed by Manning and Rosen 
[29-33]. With an approximation to centrifugal term, we solve the D-dimensional Schrodinger 
equation to its bound states energies and wavef unctions. This potential is defined as [29-33] 
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where A and a are two-dimensionless parameters [27,28] but the screening parameter b has 
dimension of length which has a potential range 1/6. The potential (1) may be further put 
in the following simple form 

Vir) = -— % , C = A, D = -A- a(a-l), (2) 

which is usually used for the description of diatomic molecular vibrations [38,39]. It is also 
used in several branches of physics for their bound states and scattering properties. The 
potential in (1) remains invariant by mapping a ^ 1 — a and has a relative minimum 
value V{ro) — — 4^^1,2^(0-1) at Tq = 61n 1 + for a > to be obtained from the first 

= 0. The second derivative which determines the force constants at r = ro 

- r=ro 

is given by 



derivative ^ 

ar 



_ A^[A + 2a{a-l)Y 

r=ro ^ ' 

The contents of this paper are as follows: In Section II we outline the Nikiforov-Uvarov (NU) 
method. In Section III, we derive I ^ bound state eigensolutions of the D-dimensional 
Schrodinger equation for the Manning-Rosen potential by this method. In Section IV, we 
present our numerical calculations in 2D and 4/} systems for various quantum numbers n 
and /. Section V, is devoted to for two special cases, namely, / = and the Hulthen potential. 
The concluding remarks are given in Section VI. 



II. THE NIKIFOROV-UVAROV METHOD 

The NU method is based on solving the second-order linear differential equation by 
reducing it to a generalized equation of hypergeometric type [37]. In this method after 
employing an appropriate coordinate transformation z = z{r), the Schrodinger equation can 
be written in the following form: 

^^^^^ + ^^"^^^ + ^^"^^^ = °' 

fpn{z) ^ (l)n{z)yn{z): (4) 
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where a{z) and a{z) are the polynomials with at most of second-degree, and f{s) is a first- 
degree polynomial. The special orthogonal polynomials [37] reduce Eq. (4) to a simple 
equation of the following hypergeometric type: 

a{z)y:{z) + T{z)y'^{z) + Xy^iz) = 0, (5) 

where 

t{z) = f{z) + 2TT{z), t'{z) < 0, 

(6) 

and A is a constant given in the form 

A = A„ = -nr'iz) - ~ ^K " {z) , n = 0, 1, 2, ... 

X^k + 7r'{z). (7) 

It is worthwhile to note that A or A„ are obtained from a particular solution of the form 
y{z) — yn{z) which is a polynomial of degree n. Further, yniz) is the hypergeometric-type 
function whose polynomial solutions are given by Rodrigues relation 

where Bn is the normalization constant and the weight function p{z) must satisfy the con- 
dition [37] 

■^w{z) = ^^w{z), w{z) = a{z)p{z). (9) 
az (^[z) 

In order to determine the weight function given in Eq. (9), we must obtain the following 
polynomial: 



7r{z) = ± 



^^'^ J -Hz) + ka(z). (10) 



In principle, the expression under the square root sign in Eq. (10) can be arranged as the 
square of a polynomial. This is possible only if its discriminant is zero. In this case, an 
equation for k is obtained. After solving this equation, the obtained values of k are included 
in the NU method and here there is a relationship between A and k given in Eq. (7). 



III. BOUND-STATE SOLUTIONS FOR ARBITRARY L-STATE 



In this section, we follow closely the approach of Ref. [36]. We begin by considering the 
SE, in arbitrary dimension D, as [40-42] 

vup[K^-nr)]}W^.c?(x) = o, 



.2 _ ^ {D-l) d 



d 



sin^-^^^_i dOD-i 



sm 



,D-2 . 



d 



^D-l 



7"2 



09 



D-1 



D-l 



.(0 



(11) 



where the potential V{r) is taken as the Manning- Rosen form in (1). In addition, /i and 
Eni denote the reduced mass and energy of two interacting particles, respectively, x is a 
D-dimensional position vector with the hyperspherical Cartesian components xi,X2, - ■ ■ ,xd 
given as follows [43-47] : 

Xi — r cos 9i sin 02 - ■ ■ sin 9d-i, 
X2 = r sin 9i sin 92 - ■ ■ sin 6'd_i, 
Xs — r cos 92 sin ^3 • • • sin 9d-i, 



Xj — r cos 9j-i sin 9j • ■ • sin 9d-i-i 3 < j < — 1, 



xd-1 = rcos6'i3_2sin6'£)_i. 
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D 

xd = rcosOn-i-, (^2) 

for D = 2, 3, • • • . We have xi = rcosip, X2 = rsimp for D = 2 and xi = r cosy? sin ^, 
X2 — r simp sin 9, x^ — rcosO for D — 3. The Laplace operator V|, is defined by [48] 

V^ = E^- (13) 

j=i ^'^i 

The volume element of the configuration space is given by 

D D-l 

Y[ dxj = r^-^drdn, dn = Y[ {sin9jy-^d9j, (14) 

i=i i=i 

where r e [0, oo), 9i e [0, 27r] and 9j e [0, tt], j e [2, D — 1]. Equation (13) permits a solution 
via separation of variables, if one writes generahzed spherical harmonics y}^l.ij^_^{^) [49]: 

= Y{h, k,---, Id-2, \m\ for D ^ 2, 

y^-in-.i^ = ^1' ^2, • • • , ^D-i) = n' H{9,), (15) 
as a simultaneous eigenfunction of : 

i^?^&._.(x)=m^y«,^_^(x), 

= ^.(^. + J - J e [1,^ - 1], 

l^O,l,---,lk^O,l,---,lk+i, ke [2,L'-2], 
h — —h, —h + 1, ■ ■ ■ ) ^2 ~ 1, ^2, 

Ll-,Yl\j5.) = l{l + D- 2)F«,^_,(x). (16) 

The unit vector along x is usually denoted by x = x/r. Additionally, the angular momentum 
operators L| are defined as [43-47,49] : 
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2 ' 



a<b=2 



1 d 



sin'=-i Ok d9k 



Lab = -i 



d__ d 



(17) 



dXi dxa 

The substitution of Eqs. (13) and (15)- (17) into Eq. (11) allows us to obtain, via the method 
of separation of variables, the following equation: 



(D-1) d 1 1 
+ - + 



dr^ r dr sin OD-idOn-i 



9 (-0-2^ ^ 



7'2 



860-1 J sin^ 9d-i 



2// _ a{a - 1) e-^"^/^ A e""/'' _ l{l + D - 2) 
+ J^Er^i (1 _ ^-r/bf + P 1 - e-l^ 



9 ( . D-2n 9 

sm^ Od-1 



2 n ' 



-L 



1 



^-2 ~ „;„D-3 



SIu''-Hd-2 900-2 



9 f ■ D-3n 9 

'sm^ -^60-2 



7'2 

^D-3 



d9o-2 ) sin^ do-2 ' 



7-2 _ 

'~ der 



(18) 
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where L|, k E [1,D — 1] are the angular operators. Thus, with the aid of Eqs. (15) and 
(16), the last wave equation can be easily separated into the following radial and angular 
parts as [50,51]:^ 

{^2 + ^pni-V,ff{r)]9ir) = 0, 



^e//(r) = ^ 



a 



(1 _ e-^/'') 1 - e-"^/^ 



+ 



+ -2)2 - 1 



(19) 



1 



sin^-2 Od-i dOo-i 



^sin^-2 Od-i^J—] +l{l + D-2) 



A 



D-2 



sin Ob-i 



H{9d-i) = 0, 



(20) 



sin^-^ 9j dOj \ ' dOi 



A._ Vl 
' sin^e, 



H{e,) = 0, je[2,L>-2], 



(21) 



H{9i) =0, Ai = l'i = 



(22) 



where Ap = lp{lp + p — \), p E [1,Z5 — 1] are separation constants and g{r-) is defined in 
Eq. (11). The solution in (22) is periodic and must satisfy the periodic boundary condition 
H{6i = </') = H{6i = </? + 27r) from which we obtain [42] 

1 



H{e,) 



exp(±iZiei), Zi = 0,1,2, 



(23) 



^The Schrodinger equation in the presence of this potential is separable to {D — \)- angular 
equations for the angular parameters (^i = (p,62, ■■■,6d-i = 0) and one radial equation for the 
radial parameter r with calculated separation constants Ap where p G [1, D — 1]. 



Further, Eqs. (20) and (21) representing the angular wave equation become 

cos6'£i_i dHiOn-i) 



+ (D-2) 



sin 6 



D-l i^^D-l 



+ 



l{l + D-2) 



A 



D-2 



D-l 



H{9d-i) = 0, 



d0] 



COS 9 J dH {9 j 



sin 9j d9j 



+ Ui 



sin^ 9i 



H{9j) = 0, 



(24) 



(25) 



with j e [2,D — 2], D > 3 and Ap which is well-known in three-dimensional space [48].^ 
Hence, Eqs. (24) and (25) are to be solved in the following subsection. 

A. The solutions of the I?-dimensional angular equations 

In order to apply NU method, we introduce a new variable s — cos9j. Hence, Eq. (25) 
is then rearranged in the form of the universal associated-Legendre differential equation 

d^H{s) js dH{s) Aj- - Aj_i - Ajs2 



+ 



-H{s) = 0, 



(26) 



ds^ ds (1-S2)2 

where j e [2, D — 2], D > 3. By comparing Eqs. (26) and (4), the corresponding polynomials 
are obtained 



r(s) = —js, a{s) = 1 — s^, a{s) = —Ajs^ + Aj — Aj_i. 

Inserting the above expressions into Eq. (10) and taking a'{s) - 
following function: 



(27) 



-2s, one obtains the 



\ 



J -2 



+ Ai-k 



s2 + A;- A, + A,_i. 



(28) 



Following the method, the polynomial 7r(s) is found to have the following four possible 
values: 



^Ad-2 = and Id-2 = m iov D = 3. 
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7r(s) = < 



(29) 



+ Aj_i) s for ki = Aj - Aj.i, 
- Aj-i) s for ki = Aj - Aj_i, 

+ A,_i for k2^Aj+[i^)\ 

- A,_i for k2 = Aj + [2^) ' , 

where Ap — Ip + {p — l)/2, with p — j — 1, j and j e [2, D — 2], D > 3 . Imposing the 
condition r'(s) < in Eq. (6), one selects the foUowing physically valid solutions: 



ki — Aj — Aj_i and 7r(s) = 



J -2 



Aj_i s, 



(30) 



which yields from Eq. (6) that 



t(s) = -2(1 + Aj_i)s. 



(31) 



Making use from Eq. (7), the following expressions for A are obtained as follows: 



A = = 2nj{l + Aj_i) + nj{nj - 1), 



(32) 



X = Aj - Aj_i - Aj_i 
Upon comparing Eqs. (32) and (33), we obtain 



J -2 



(33) 



(34) 



In addition, using Eqs. (6) and (8)-(9), we obtain the following useful parts of the wave- 
functions: 



(35) 



where j e [2, D — 2], D > 3. Besides, substituting the weight function p{s) given in (35) into 
Eq. (8), we obtain 

,nj+Aj_i 



where is the normaliation factor. Finally the angular wavefunction is 
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(36) 



Hn,{e^) = A^n, (sin^,)'-^ Pif-^'^-^)(cos^,), J e [2, L> - 2], L> > 3 



(37) 



with rij given in Eq. (34) becomes 

Uj = Ij - Ij-i, j e[2,D -2\,D > 3. (38) 

Likewise, using s — cos9d-i, we can rewrite Eq. (24) in the associated Legendre form 

d^H{s) {D-l)sdH{s) u{l - s^) - Ad-2 



1 — ds 



+ 



'-H{s) = 0, 



(39) 



(40) 



It's worth to note that, Eq. (39) has been recently solved in 3D by the NU method in 
[40,41,51]. However, our aim is to solve it in £)-dimensions. Hence, comparing Eqs. (39) 
and (4), the corresponding polynomials are obtained 



f{s) = -{D-l)s, cr(s) = l-s^ a{s) = -us'^ + u - A 



D-2- 



(41) 



Inserting the above expressions into Eq. (10) and taking a'{s) = —2s, one obtains: 



D-3 



+ iy-k 



D-2- 



(42) 



Following the method, the polynomial 7r(s) is found to have the following four possible 
values: 



7r(s) = < 



^ + Ad-2) s for /ci = 1/ - Ad-2, 

- Ad-2) s for /ci = i/ - Ad-2, 

i£^s + An-2 iork2^u+{^)\ 

i^s -An-2 iork, = iy+ {^)\ 



(43) 



where A0-2 — Id-2 + ^2^- Imposing the condition t'(s) < in Eq. (6), one selects the 
following physically valid solutions: 



ki = u — Ad-2 and 7r(s) = —Id-2S, p = l{l + D — 2), 



(44) 
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giving 

t(s) = -2(1 + Ao_2)s. (45) 
Making use from Eq. (7), we obtain 

A = A„o_^ = 2nD-i{l + Ad-2) + nD-i{nD-i - 1), (46) 

X = l{l + D-2)-lD-2{iD-2 + D-2). (47) 

We compare Eqs. (46) and (47), the angular momentum / values are obtained from 

/ = hd-i + Id-2, (48) 

which can be easily reduced to the well-known solution 

I = n + m, (49) 

in 3D [51]. Using Eqs (6) and (8)-(9), we obtain the following useful parts of the wavefunc- 
tions: 

0(.) = (i - .^y^-^/^ p(.) = (i - (50) 

Besides, the Rodrigues relation (8) gives 

yn,.,{s)^B^^_^(l-s') -^—-(l-s') , (51) 

where -B„^_j is the normaliation factor. Finally the angular wavefunctions are 

Hn,_,{eD-i) = Nn,_, (sine^_i)'--^ Pi^--'^--)(cosez?_i), (52) 
where tid-i is given by Eq. (48). 
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B. The solutions of the D-dimensional radial equation 



Since the radial part of the D-dimensional Schrodinger equation with above Manning- 
Rosen effective potential has no analytical solution for I ^ states, an approximation to the 
centrifugal term has to be made. The good approximation for in the centrifugal barrier 
is taken as [18,33]^ 



-r/b 



2 ' 



(53) 



in a short potential range. To solve it by the present method, we need to recast Eq. (19) 
with the aid of Eq. (53), into the form of Eq. (4) changing the variables r — > z through the 
mapping function r = f{z) and energy transformation given by 



z = e-'l\ e 



72 ' ^nl < U, 



(54) 



to obtain the following hypergeometric equation: 



d'g{{z) ^{1-z) dg{z) 



dz^ 



z{l — z) dz 



1 



[z{l - z)] 



4^ + - (D + 2/ - 2)2 + 1 z- A + e^ + a{a - 1) z^\ g{z) = 



(55) 



We notice that for bound state (real) solutions, the last equation requires that 



(56) 



0, when r — > oo, 

1, when r ^ 0, 

and thus the finite radial wavefunctions Rni{z) — > 0. To apply the NU method, it is necessary 
to compare Eq. (15) with Eq. (4). Subsequently, the following value for the parameters in 
Eq. (4) are obtained as 



1 e-2'-/*' 



The series approximation to the expression vi-r^ 

( 1— 



7 « p — ^, it includes a Coulomb term. 
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t{z) — 1 — z, a{z) — z — 



a{z) = - A + e^ + a{a - 1) 



AA + 8e^ -{D + 21- 2f+ 



z — e . 



(57) 



If one inserts these values of parameters into Eq. (10), with (j'{z) = 1 — 2z, the following 
linear function is achieved 

7r(^) = -J 



±^^{1 + 4 [A + £2 + q;(q; - 1)] - A;} ^2 + {Ak - [4A + Ss^ _ (£) + 2Z - 2)2 + 1]} z + 4£2. 

(58) 

According to this method the expression in the square root has to be set equal to zero, that is, 
A = {1 + 4 [A + £2 + q;(q; - 1)] - A;} ^2 + {4A; - [4A + Ss^ _ (£) + 2Z - 2)2 + 1]} z + 4£2 = 0. 
Thus the constant k can be obtained as 



^ ^ ^ _ (D + 2[_^)2 — 1^^^^ a = <^(1 - 2a;)2 + (L» + 2Z - 2)2 - 1. 

In view of that, we can find four possible functions for 77(2;) as 

for k^A- (^+2^-2)^-1 + as, 
for ^ A - (^+2^-2)^-1 _ 



(59) 



(60) 



We must select 

^ (L> + 2Z- 2)2-1 ^ / a\ 

k = A ^ as, 7r(^) = "2 ^ ^ ~ V ^ 2J ^' 

in order to obtain the polynomial, t{z) = f{z) + 2tt{z) having negative derivative as 
t{z) = l + 2e-{2 + 2e + a)z, t'{z) = -(2 + 2e + a). 



(61) 



(62) 



We can also write the values of A = A; + tt'{z) and A„ = —nr'^z) — ^\ "{z), n = 0, 1, 2, ... 



as 



X = A-^ ^-^ (l + «) 



(63) 
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A„ = n(l + n + a + 2s), n = 0, 1, 2, ... (64) 

respectively. Additionally, using the definition of A = A„ and solving the resulting equation 
for £, allows one to obtain 

4(n + 1)2 + (D + 2/ - 2)2 - 1 + 4(2n + l)r; - 4yl -1 + a 



from which we obtain the discrete energy levels 



2 



32/i62 



4(n + 1)2 + (D + 2/ - 2)2 + 4(2n + l)'q-AA-l 



2 



, < n, Z < oo (66) 



2(n + 1 + 77) 

where n denotes the radial quantum number. It is found that A remains invariant by 
mapping a ^ 1 — a, so do the bound state energies Eni. An important quantity of interest 
for the Manning-Rosen potential is the critical coupling constant which is that value of 
A for which the binding energy of the level in question becomes zero. Using Eq. (26), in 
atomic units — ^ — ^ — e = l, 

.le = (n + l + 77)^-77(77 + l) + ^^±|^-^. (67) 

Let us now find the corresponding radial part of the wave function. Using (t{z) and 'k{z) 
in Eqs (57) and (61), we obtain 

(t>{z)^z'{l-zi''+''^/\ (68) 
p{z)=z^'{l-zf'^+\ (69) 

y^i{z) = C„^-2^(l - z)-^'"^'^^ [z^+'^l - zT^'^^'] . (70) 

The functions yni{z) are, up to a numerical factor, are in the form of Jacobi polynomials, 
i.e., yni{,z) ~ p^2£,2r;+i)(^]^ _ 2^;)^ vafid physically in the interval (0<r<oo^0<2;<l) 
[52] . Therefore, the radial part of the wave functions can be found by substituting Eqs. (68) 
and (70) into Rni{z) = (f){z)yni{z) as 

9ni{z) = A^„;^^(l - zy^^Pi'^''^^'\l - 2z), (71) 
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where e and A are given in Eq. (65) and Nni is a normalization constant. This equation 
satisfies the requirements; Rni{z) = as z = (r ^ oo) and Rni{z) = as 2; = 1 (r = 0). 
Therefore, the wave functions, gni{z) in Eq. (71) is valid physically in the closed interval 

2;e[0,l]orre(0,oo). Further, the wave functions satisfy the normalization condition 

00 1 

/ \9ni(r)fdr ^l = bjz-^ \gni(z)fdz, (72) 


where Nni can be determined via 

1 

1 = bNli J z^^-\\ - z)'^+^ [Pi^'''^+'\l - 2z)]'dz. (73) 


Following Ref. ]36}, we find the normalization constant 

Nni = -1^, (74) 
^s(n) 



where 



_ r(n + 2^ + 2)r(n + 2e + l)^ 

'^'')-^^ ^) r{n + 2s + 2r] + 2) 



^ " (-l)P+T(n + 2e + r -p+l){p + 2ri + 2) 

^ pj:^oP\r\{n-py.{n-ry.T{n + 2e-p+l)r{2e + r + l){n + 2e + r + 2r] + 2)' ^ ^ 

Therefore, we may express the normalized total wave functions as 

V'(X) = Ar„;r-(^-3)/2g-sr-/6^^ _ ^-r/bY+r,pi2e,2n+l)^^ _ ^^-r/b^ exp(±iZi^i) 

D-2 

xNn,_, sin(^^_i)'--Pi^^--^--)(cos^,,_i). n {sme.y^-' Pif-^'^-^Hcos^,-) (76) 

To show the accuracy of our results, we calculate the energy eigenvalues for various n and 
/ quantum numbers with three different values of the parameters a. As seen n Table 1, the 
energy eigenvalues for different quantum numbers are obtained numerically for D = 2 and 
D — A cases. 

/^From Eq. (66), we have seen that two interdimensional states are degenerate whenever 

[53] 
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{n,l,D) ^ {n,l±l,D^2). (77) 

Thus, a knowledge of ior D — 2 and D — 3 provides the information necessary to find 
E^'' for other higher dimensions. 

For example, Eq^I = Eq^I = E^^l = -^o^i- This is the same transformational invariance 
described for bound states of free atoms and molecules [54,55] and demonstrates the existence 
of interdimensional degeneracies among states of the confined Manning-Rosen potential. 

As an example of incidental degeneracy. Table 1 presents the results of the confined 2D 
and 4D Manning- Rosen energies at several radii of confinement for various n and / states. 



IV. DISCUSSIONS 

In this work, we have utilized NU method to solve the £)-dimensional SE for the Manning- 
Rosen model potential with the angular momentum I ^ states. We have derived the binding 
energy spectra in Eq. (66) and their corresponding wave functions in Eq. (71). 

Let us study special cases. We have shown that for a — (1), the present solution 
reduces to the one of the Hulthen potential [16,18,19]: 

yW(^) = -^0^4^' ^0 = Ze% 6 = b-' (78) 

where Ze^ is the strength and 5 is the screening parameter and b is the range of potential. 
If the potential is used for atoms, the Z is identified with the atomic number. This can be 
achieved by setting rj — ^{D + 21 — 3), hence, the energy for Z 7^ states 

[4:A-(2n + D + 2l-lYfn^ 

and for s-wave (/ = 0) states 

Essentially, these results coincide with those obtained by the Feynman integral method [31] 
and the standard way [32,33], respectively. Furthermore, if taking b — 1/5 and identifying 
as Ze^(5, we are able to obtain 

17 



2^2 



■(2n + L'H-2/- 1) 



(2n + D + 2Z-l) %Ze'^ix 
which coincides with those of Refs. [16,18]. With natural units = ji = Z 
have 

^ [ 1 (2n + £> + 2Z-l) 1^ 

[(2n + D + 2Z-l) 8 

which coincides with Refs. [16,33]. 



(81) 

1, we 

(82) 



The corresponding radial wave functions are expressed as 



Rnl{r) = N^ie-''\l - e-^'-)(^+2^-l)/2p^2e,M-2)(^ _ ^^-5ry-(D-m ^ 



(83) 



where 



2/iZe2 



n^5 



n^5 [(2n + D + 2Z-l) SZe^n 



{2n + D + 2l- 1) 



, < n, / < oo. 



(84) 



which coincides for the ground state with that given in Eq. (6) by Goniil et al [18]. In 
addition, for 5r <^ 1 (i.e., r/h <S 1), the Hulthen potential turns to become a Coulomb 
potential: V{r) — —Ze^/r with energy levels and wavef unctions: 

4£o 



E.. 



nl 



(2n + D + 2/ - 1)2 



, n = 0,l,2,.. 



2*2 



_ z'h 



11& 



(85) 



where Sq — 13.6 eV and ao is Bohr radius for the Hydrogen atom. The wave functions are 



Rni = Nni exp 



8/xZe2 



r {2n + D + 2l- 1) 



^{D+2l-l)/2p{.f 



,D+2l-2 



(1 + 2Sr) 



which coincide with Refs. [3,16,22]. 



V. COCLUDING REMARKS 



In this work, we have extended the approximate solutions of the /-wave Schrodinger 
equation with the Manning-Rosen potential to D-dimensions . The special cases for a = 0, 1 
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are discussed. The results are found to be in good agreement with those obtained by other 
methods in 3D for short potential range, small a and / [36]. These numerical solutions have 
also extended to various dimensional space, D = 2 and D = 4 systems. We have also 
studied two special cases for Z = 0, Z 7^ and Hulthen potential. The results we have ended 
up show that the NU method constitute a reliable alternative way in solving the exponential 
potentials. 
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TABLES 

TABLE I. Eigenvalues for 2p, 3p, 3d, 4p, 4(i, 4/, 5p, 5d, 5/, Sgi, 6p, 6d, 6/ and 6^ states in atomic 
units {Ti = II = 1) and for a = 0.75 and a = 1.5, A = 26. 

i:) = 2^ i:) = 4^ 

states 1/6 a = 0.75 a = 0, 1 a = 1.5 a =0.75 a =0,1 a = 1.5 

2p 0.025 -0.241087728 -0.209898003 -0.140949065 -0.070734690 -0.067988281 -0.058898861 
0.050 -0.227946676 -0.197925347 -0.131737328 -0.059344084 -0.056953125 -0.049054156 
0.075 -0.215173874 -0.186304253 -0.122836866 -0.048952839 -0.046894531 -0.040109106 
0.100 -0.202769319 -0.175034722 -0.114247678 -0.039560954 -0.037812500 -0.032063712 

3p 0.025 -0.074279113 -0.067988281 -0.051933432 -0.030209821 -0.029273358 -0.026068346 
0.050 -0.062813564 -0.056953125 -0.042142549 -0.020395577 -0.019644452 -0.017092049 
0.075 -0.052308602 -0.046894531 -0.033373420 -0.012502916 -0.011929608 -0.010003237 
0.100 -0.042764227 -0.037812500 -0.025626042 -0.006531840 -0.006128827 -0.004801908 

3d 0.025 -0.070734690 -0.067988281 -0.058898861 -0.029833656 -0.029273358 -0.027228277 
0.050 -0.059344084 -0.056953125 -0.049054156 -0.020047209 -0.019644452 -0.018176769 
0.075 -0.048952839 -0.046894531 -0.040109106 -0.012199670 -0.011929608 -0.010947973 

Ap 0.025 -0.031448122 -0.029273358 -0.023381941 -0.014180352 -0.013773389 -0.012357598 
0.050 -0.021545731 -0.019644452 -0.014606136 -0.006296995 -0.006019483 -0.005072360 
0.075 -0.013510134 -0.011929608 -0.007885467 -0.001570215 -0.001429639 -0.000978205 

Ad 0.025 -0.030209821 -0.029273358 -0.026068346 -0.014011823 -0.013773389 -0.012892982 
0.050 -0.020395577 -0.019644452 -0.017092049 -0.006162813 -0.006019483 -0.005494347 
0.075 -0.012502916 -0.011929608 -0.010003237 -0.001492711 -0.001429639 -0.001204122 

4/ 0.025 -0.029833656 -0.029273358 -0.027228277 -0.013929374 -0.013773389 -0.013182139 
0.050 -0.020047209 -0.019644452 -0.018176769 -0.006097355 -0.006019483 -0.005724889 
0.075 -0.012199670 -0.011929608 -0.010947973 -0.001455297 -0.001429639 -0.001333163 

5p 0.025 -0.014732070 -0.013773389 -0.011100961 -0.007127957 -0.006916484 -0.006175251 

5c? 0.025 -0.014180352 -0.013773389 -0.012357598 -0.006506751 -0.006392207 -0.005967020 
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5/ 0.025 

5g 0.025 

6p 0.025 

6d 0.025 

6/ 0.025 

6g 0.025 



-0.014011823 
-0.013929374 
-0.006866319 
-0.005435481 
-0.006506751 
-0.006465489 



-0.013773389 
-0.013773389 
-0.006392207 
-0.006392207 
-0.006392207 
-0.006392207 



-0.012892982 
-0.013182139 
-0.005056211 
-0.005695750 
-0.005967020 
-0.006113207 



-0.006465489 
-0.006440958 
-0.002734814 
-0.002691847 
-0.002670817 
-0.002658317 



-0.006392207 
-0.006392207 
-0.002635101 
-0.002635101 
-0.002635101 
-0.002635101 



-0.006113207 
-0.006204004 
-0.002286461 
-0.002424502 
-0.002499036 
-0.002545374 



^Two-dimensional Schrodinger equation. 
'^Four-dimensional Schrodinger equation. 
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